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Consider  a multicoirimodity  transhipment  problem  where  the  prices 

at  each  location  are  an  affine  function  of  the  supplies  and  demands 

i 

at  that  location  and  the  shipping  costs  on  a link  are  an  affine 

function  of  the  quantities  shipped  on  that  link.  A system  of  prices, 

supplies,  demands,  and  shipments  is  defined  to  be  an  equilibrium,  if 

I there  is  a balance  in  the  shipments,  supplies,  and  demands  of  goods 

I 

at  each  location,  if  local  prices  do  not  exceed  the  cost  of  Importing, 

‘ , and  if  shipments  are  price  efficient.  Lemke's  algorithm  is  used  to 

compute  an  equilibrium  in  a finite  number  of  steps. 

This  paper  extends  the  works  of  Takayama  and  Judge  [4]  who 
utilized  quadratic  programming.  Our  approach  and  results  are  more  direct 
and  more  general;  we  solve  the  equilibrium  conditions  directly  without 
j passage  to  an  optimization  problem.  In  [3]  an  even  more  general  equlli- 

^ brium  problem  is  formulated  as  a complementarity  problem,  but  in  a dif- 

ferent manner  than  that  used  here  and  without  the  results  obtained  here. 

Our  problem  is  conveniently  represented  with  a directed  graph 
(..f',  j/’)  with  a finite  number  of  nodes  .A'*  {1,  ...,  n}  and  a finite 
number  of  links  i?*  {1,  ...,  i).  Such  a directed  graph  with  n - 4 


/ 

Figure  1 


Each  node  1 In  represents  a producer/consumer  at  a 
specific  spatial/temporal  location.  Each  link  s in  ^ represents 
a specific  transport  facility  for  transfering  comnodities  between 
nodes,  that  is,  locations,  and  each  link  s is  oriented  to  coincide  with 
the  direction  of  a possible  transfer  activity.  For  example,  if  link  s 
has  head  J and  tall  i,  that  is. 


then  the  transfer  of  goods  along  a is  from  node  1 to  node  J . 


We  permit  multiple  links  between  the  same  pair  of  nodes, e.g.. 


I 


! 

I 

I 

I 

I 

I 


I 


1 

1,  2 A' 

1,  3 e 

but  we  do  not  permit  loops,  e.g., 

1 e.yf' 


(Technically  speaking,  ^ indexes  a finite  subset  of  {(i,  j , z)  ; 

i,  j e ; i 5^  j ; z = 1,  2,  . ..}.) 

We  are  concerned  with  the  supply,  demand,  and  transhipment  of 

m goods  g = 1,  ...,  m in  the  network  A good  g could  be 

a raw  material,  an  intermediate  product,  or  a finished  product. 

Let  the  variable  f * (f  ,,...,  f ) represent  the  quantities 
s s X sm 

of  the  various  goods  shipped  along  link  s and  let  the  per  unit  shipping 

cost  of  the  various  goods  moved  along  link  s be  an  affine  function 

f + c . Hence,  f • (C  f + c ) is  the  total  shipping  cost  on 
s s s s s s s 

link  8. 


3 


Let  the  variable 


Pj  “ represent  the  prices 


of  the  various  goods  at  node  1.  Let  the  variable  h.  • (h.,,  h.  ) 

1 11  in 

represent  the  net  exports  from  node  1 of  the  various  goods;  that  Is, 
h^g  Is  positive  when  node  1 produces  more  of  good  g than  It  consumes, 
and  Is  negative  when  more  Is  consumed  than  produced.  We  assume  that  the 
prices  p^  and  exports  h^^  at  node  1 are  related  by  an  affine  function 
Pi  - A^h^  + a^. 

Our  network  Is  completely  specified  by  the  graph  (./f,  y) , the 
(A^la^)  for  1 in  .4',  and  the  s in  .5?.  For  the 

graph  of  Figure  1 we  could  have,  for  example,  the  data  of  Figure  2 
corresponding  to  two  goods;  In  Section  A we  shall  solve  this  example. 


Figure  2 
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To  define  an  equilibrium  In  a precise  manner  we  need  some 
additional  notation.  For  a link  s,  let  as  be  the  node  at  the  tall 
of  s,  and  let  sa  be  the  node  at  the  head  of  s.  For  example,  with 
respect  to  Figure  1,  we  have  al  = 1 and  5n  * 2.  For  a node  1,  let 
be  the  set  of  links  s entering  1,  and  i-»^  be  the  set  of  links 
leaving  1.  Hence,  link  s is  in  !-»  or  >1  if  and  only  if  ns  = i 
or  sa  =1,  respectively.  With  respect  to  Figure  1,  we  have  !->■  = {1} 
and  *-2  = {1,5}. 

A system  of  prices  p = (p, , . . . , p ) , exports  h = (h, , . . . , h ) , 

in  in 

and  shipments  f = •••.  definition  an  equilibrium  if  the 

following  five  conditions  hold. 


(a)  fg  > 0 


s e 


(W  I - I 


1 -*■  ->1 


1 


(1) 


(c)  p^  » A^h^  + a^ 


i 


(d)  p +Cf+c>p 

ns  3 s s — sa 


s G </’ 


(e)  f*(p  +Cf  +c  -p  )=0 
s ns  s s s sa 


Condition  (la)  requires  that  the  shipments  or  flows  be  nonnegative, 
however,  note  that  the  prices  and  exports,  may  be  positive,  negative. 
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or  zero.  Condition  (lb)  represents  the  conservation  of  goods  at 

node  1.  Condition  (Ic)  expresses  the  fact  that  node  1 produces 

and/or  consumes  according  to  the  prices.  Condition  (Id)  Is  a price 

stability  condition  requiring  that  the  local  price  must  not  exceed 

a neighboring  price  plus  transportation  costs.  Finally, condition  (le) 

requires  that  goods  only  be  shipped  (In  positive  amounts)  on  price 

efficient  links,  that  Is  , If  f Is  positive,  then 

sg 


(p  + C f + c ) 
as  s s s g 


'Psn>g 


the  price  of  g at  as  plus  transportation  costs  along  s to  sa 
must  equal  the  price  of  g at  sa. 

In  the  following  sections  we  prove  and  Illustrate  the  following 
theorems.  Theorem  1 provides  conditions  under  which  Lemke's  algorithm 
will  find  an  equilibrium  or  show  that  none  exists.  Theorems  2 and  3 
give  additional  conditions  under  which  Lemke's  algorithm  will  always 
find  an  equilibrium. 


Theorem  1.  If  each  Is  positive  semldeflnlte  and  each  C^  Is 

copositive  plus,  then  Lemke's  algorithm  generates  an  equilibrium  or 
demonstrates  that  no  equilibrium  exists.  13 


Theorem  2.  If  each  Is  positive  semldeflnlte  and  each  C^  ‘Is 

strictly  copositive,  then  Lemke's  algorithm  generates  an  equilibrium.  El 
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u 

Theorem  3.  If  each  is  positive  definite  and  If  for  each  link 

s,  C is  copositive  pins  and  t > 0,  i -C  f = 0,  c *1  <0  has 

s S“  SSS  ss 

no  solution,  then  Lemke's  algorithm  generates  an  equilibrium.  0 


i 

} 

i 

I 

I 

) 

1 


- 1 


Note  that  sn-^  Is  the  set  of  all  links  whose  tails  are  the  head 
of  s,  etc.  Hence,  If  we  solve  (3)  and  use  (lb)  and  (Ic)  to  compute 
hj  and  we  have  solved  (1)  and  have  an  equilibrium;  conversely, 

any  solution  to  (i)  yields  a solution  to  (3). 

Equation  (3a)  can  be  rewritten  as 


w =(A  +A  +C)f 
s ns  sn  s s 


+ (A  +A)yf  -(A  + ) I 

ns  sn  ^ t ns  sn  J t 


(4) 


+ A If, 

ns  ^ t 
Y 


A I f 

ns  4 t 


+ A If 
sn  , t 


-A  I f,  + 
sn  ^ t 


where 


a = ((ns*-)  (-(-sn))  'v  s 6 * (sn-v)  H (-ms) 


Y = (ns-^)  'V  (mn)  6 “ (ms)  (sn-*-) 


^ “ (mn)  'V  (ns->-)  n ••  (sn-*-)  'v  (ms) 
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Note  that  the  index  sets  a through  n are  pairwise  disjoint;  the 
following  schema  illustrates  the  partition: 


Now  let  us  cast  the  system  (3)  into  the  linear  complementarity 
problem.  Define  an  Jim  x S-m  partitioned  matrix  M by  M « with 

s,  t “ 1,  H and  is  an  m x m matrix  defined  by 


Define  an  tm  vector  q to  be  (c^,  . Using  (4)  it  is  now 

evident  that  solving  (2)  is  equivalent  to  solving  (3). 

As  an  example,  consider  the  multicommodity  network  of  Figure  1, 
we  have  displayed  the  matrix  (M|q)  in  Figure  3. 
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3.  SOLUTION 


We  now  prove  the  three  theorems  stated  In  Section  1.  A square 

matrix  N Is  said  to  be  positive  definite  If  xNx  Is  positive  for 

all  nonzero  x.  N is  said  to  be  positive  semldeflnlte  If  xNx 

Is  nonnegative  for  all  x.  N Is  said  to  be  strictly  copositive  If 

xMx  is  positive  for  all  nonzero  nonnegative  x.  N Is  said  to  be 

copositive  plus.  If  xNx  Is  nonnegative  for  all  nonnegative,  x and 

T 

if  xNx  equals  zero  for  nonnegative  x Implies  (N  + N)x  equals 
zero.  Incidently,  for  clarity,  we  remark  that  the  and  M 

are  not  assumed  to  be  symmetric. 

For  any  system  of  shipments  f and  exports  h satisfying  the 
conservation  equations,  we  have  the  equality 


(5) 


f ‘Mf 


I hj^-A^h^  + I 


.A- 


•X 


•c  f 

s s 


From  our  formula  (5)  we  obtain  the  following  lemma. 


Lemma  4.  If  each  A^  Is  positive  semldeflnlte  and  each  Is  co- 

positive plus,  then  M Is  copositive  plus. 

Proof . If  each  is  copositive  plus,  then  so  Is  the  matrix 

C * dlag(C^,  ...,  . Also,  clearly  the  matrix  M > C la  positive 

semldeflnlte,  and,  therefore,  copositive  plus.  Finally,  the  sum  of 
two  copositive  plus  matrices  Is  copositive  plua.  13 
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Towards;  solving  the  linear  complementarity  problem  (2)  we 


V 


I 

L 


I 


i 

i 

I 

J 

t 

» 

< 

i 


i 

i 

4 

i 


apply  I.emke's  algorltbm  to  t!ie  system 


(6) 


Iw  - Mf  - ez  * q 

wO  f>0  wf“0  z>0 


Briefly,  Lemke's  algorithm  proceeds  by  generating  a path  of  solutions 
to  (6)  beginning  with  the  family  of  solutions  for  which  f « 0.  A 
complete  description  of  Lemke's  algorithm  can  be  found  in  Lemke  121, 
or  [1]. 


Theorem  5 (Lemke) . Suppose  M is  copositive  plus.  If  Lemke's 
algorithm  is  applied  to  (6)  whore  e - (e^,  e^,  . . . , e^,  I 

if  '^k  - ®k  “ ^ '^k  ^ either  a solution  to  (6)  with 

2=0  is  generated  or  a vector  f is  generated  satisfying  (7). 


(7)  fM  < 0 f-Mf  - 0 f-q  <0  f i 0 

Proof.  See  (IJ.  E 
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Theorem  1 now  follows  from  Lemma  4 and  Theorem  5;  note  chat  (2) 
cannot  have  a solution  If  (7)  does. 


Proof  of  Theorem  2.  Suppose  the  hyptheses  of  Theorem  2 are  satisfied 
but  Lemke's  algorithm  does  not  find  a solution.  From  (5)  and  (7) 

0 . f-Mf  . I 

Since  each  term  is  nonnegative,  each  term  is  zero.  Since  each 
is  strictly  coposltlve,  each  f^  equals  zero.  But  this  contradicts 
f*q  < 0.  So  Lemke's  algorithm  will  generate  a solution.  S 

Proof  of  Theorem  3.  Suppose  the  hypotheses  of  Theorem  3 are  satisfied, 
but  Lemke's  algorithm  does  not  find  a solution.  From  (5)  and  (7) 

0 - t-Mf  . I + I 

Since  each  term  is  nonnegative,  each  term  is  zero.  Each  h^  is  zerc 
since  each  is  positive  definite.  Also 

0 > y f_*(c  +a  -a  )”yc*f  +y  a, ‘h,  * ) c 

‘•8  8 ns  Sn  ‘*88  ^ii  ^ 88 

so  for  some  link  s,  c ‘f  < 0,  f *0  f - 0,  f >0  which  contradicts 

SB  BBS  S — 

the  hypotheses.  So  Lemke's  algorithm  will  generate  a solution.  El 
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4 . 

We  shall  now  Illustrate  the  preceding  results  by  solving  the 
2-comn;o<jity  network  problem  defined  by  the  graph  of  Figure  1 and 
the  data  of  Figure  2.  (M,  q)  is  shown  In  Figure  4,  see  Figure  3. 

I.emke's  algorithm  is  applied  as  described  In  Theorem  5.  In  view  of 
Theorem  2 convergence  to  an  equilibrium  Is  guaranteed;  the  flow 
generated  Is 


Hence,  (p,  h,  f)  Is  an  equilibrium  when  the  prices  p and 
exports  li  are  computed  according  to  (lb)  and  (Ic) . 


15 


REFERENCES 


1 

I 

j 


1.  EAVES,  B.  C.,  The  Linear  Complementarity  Problem",  Management 

Science,  17,  9,  May  1971,  612-634. 

2.  l.EMKE,  C.  E. , Blmatrlx  Equilibrium  Points  and  Mathematical  Pro- 

grammtng.  Management  Science,  11,  1965,  681-689. 

3.  PETERSON,  E.  1..,  The  Conical  Duality  and  Complementarity  of  Price 

and  Quantity  for  Multicommodity  Spatial  and  Temporal  Network 
Allocation  Problems,  Center  for  Math.  Studies  in  Economics 
and  Management  Science,  Disc.  Paper  #207,  Northwestern 
University,  Evanston,  111.,  March  1976. 

4.  TAKAYA.MA,  T.,  and  G.  G.  JUDGE,  Spatial  and  Temporal  Price  and 

Allocation  Models,  1971,  Nor th-Holland  Pub.  Co. 


I 


1 

I 


i 

i 


17 


UNCIJ^SSIH 


SCCUJBiTV  CL*SSir|C*TlON  O*  This  l'*TT<«n  Dmit  Inlmttdy 


REPORT  DOCUMENTATION  PAG 


HtFORF.  COMPLETINC.  KOFfM 


J MlCl^lEllT  s CATALOC  NUMdtH 


^ GOVT  AcceiiiON  NO 


TtTLt 


TYPE  or 


Computing  Economic  Equilibria  on  Affine  Networks 
with  Lcmke's  Algorithm  | 


Technical 


^INACT  on  GMANT  M<IMacn|b; 

■AFOSR-yT-lllA-  . I 


sniuth 


DAA(Q9-74-C-//3Z 


* PCArORMING  organization  name  ano  address 


Program  element  rroject.  task 

AREA  A WORK  UNIT  NUMBERS 


Department  of  Operations  Research 
Stanford  University 


II.  CONTROLLING  OrriCE  NAME  AND  ADDRESS  i 

Air  Force  Office  of  Scientific  Research  ( 
Air  Force  Systems  Command,  United  States  Air  Fore' 
Boeing  Air  Force  Base,  Washington,  D.C.  20332 


AGES 


II.  SECU 


Mathematics  Division 
U.S.  Army  Research  Office 
Box  CM,  Duke  Station 
Durham,  North  Carolina  27706 


UNCIASSIFTED 


OECLASSiriCATION/ DOWNGRADING 
SCHEDULE 


II  DISTRIBUTION  statement  (oI  thit  Ktpotl) 


Approved  for  public  release  and  sale:  distribution  is  unlimited 


t7.  OlSTmauTlON  STATCMCNT  (oi  th0  •bttemet  In  VfocA  90,  It  ^IllnfOnl  Ifm 


‘ ' rv'c  I).  ' report  are  nnt  t i be  conftrued  ns  ati 
i'  ;..;.ri.,.i.iii  of  Uie  Army  posiuou,  unless  80 
iinsignated  by  uiher  uutiiui igcii  ducunienUi. 


!•  KCr  WOffOS  (Conrinw*  on  r«v«r««  0I4O  It  noco§oofr  Idonlltf  Of  OtoeA  ntmbor) 


Complementary  pivot  theory,  Economic  equilibria,  Lemke's  algorithm 
Multicommodity  networks. 


tmr  ond  Idontitf  Of  OlotO  mimOof) 


Economic  equilibria  between  supply  and  demand  on  certain  affine  multl- 
comnodlty  networks  are  characterized  as  solutions  to  a linear  complementarity 
problem  to  which  Lemke's  algorithm  is  applied. 


